Detection of the electrostatic spin crossover effect in magnetic molecules 
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Scanning tunneling microscopy (STM) can be used to detect inelastic spin transitions in mag- 
netic nano-structures comprising only a handful of atoms. Here we demonstrate that STM can 
uniquely identify the electrostatic spin crossover effect, whereby the exchange interaction between 
two magnetic centers in a magnetic molecule changes sign as a function of an external electric field. 
The fingerprint of such effect is a large drop in the differential conductance as the bias increases. 
Crucially in the case of a magnetic dimer the spin crossover transition inverts the order between the 
ground state and the first excited state, but does not change their symmetry. This means that at 
both sides of the conductance drop associated to the spin crossover transition there are two inelas- 
tic transition between the same states. The corresponding conductance steps split identically in a 
magnetic field and provide a unique way to identify the electrostatic spin crossover. 



Probing and manipulating spins in solid state systems 
underpin the future development of spintronics and quan- 
tum information technology. In particular the interaction 
between conduction electrons and transition metal atoms 
with partially filled d-shells drives many low-temperature 
spin effects. Analyzing the conductance spectra of mag- 
netic atoms adsorbed on a metallic host and probed by a 
scanning tunneling microscope (STM) unearths various 
distinctive features, which are indicative of many-body 
scattering. These include conductance steps associated 
to spin- flip transitions and zero-bias Kondo resonances. 
Such probing method, called spin-flip inelastic electron 
tunneling spectroscopy (SP-IETS), allows one to extract 
microscopic information about the spin interaction at the 
single atom level [IHl]. 

So far SP-IETS experiments have all been well ex- 
plained by assuming that the low-energy spin Hamilto- 
nian for the atoms to probe is insensitive of the current 
or the STM voltage. Yet, recently it has been theoreti- 
cally demonstrated that an electrostatic potential of suf- 
ficient strength can alter the exchange interaction be- 
tween magnetically-coupled atoms [TOl - fTB] . In partic- 
ular, it was shown that the magnitude of the energy 
Stark shift in a molecule containing two magnetic cen- 
ters depends on whether the two spins are either parallel 
(high-spin) or antiparallel (low-spin) to each other, lead- 
ing to a quadratic dependence of the exchange interac- 
tion between the two centers on the applied bias volt- 
age [TU]. Furthermore, when an electrical dipole breaks 
the molecule inversion symmetry, the spin crossover may 
occur at experimentally achievable electric fields. This 
effect paves the way for the development of quantum 
information and spintronics devices |14) . Preliminary 
experimental evidence of such, so called electrostatic 
spin crossover effect (ESCE), has been already provided 
in breaking junctions experiments with organometallic 
molecules |15|. These experiments however are complex 
since the molecule can access both different spin and 
charging states. 

Here we investigate the possibility of using STM to de- 



tect the dependence of the exchange coupling on an elec- 
trical potential, by using our previously developed quan- 
tum mechanical approach based on the non-equilibrium 
Green's function formalism [TBHIH]- In particular, we 
calculate the bias-dependent conductance spectra of an 
ideal molecule comprising two exchange-coupled spin 1/2 
atoms. Previous attempts to investigate the ESCE have 
used a classical description of the localised spins [19) . 
where quantized excitations were neglected. In con- 
trast our fully quantum mechanical approach allows us 
to study in details the elementary excitations of the sys- 
tem and identify the fingerprint of the ESCE. We will 
show that there exists a critical voltage, Vc, where the 
conductance profile changes drastically and that this cor- 
responds to the spin crossover between a low- and a high- 
spin state. Importantly, as the spin crossover transition 
reverses the order between the ground state and the first 
excited state but not their symmetry, at both side of Vc 
there is an inelastic conductance step between spin states 
having same total spin. These split in the same way in a 
magnetic field and provide a unique way to identify the 
ESCE-driven transition. 

We consider a single-orbital tight-binding model [TBl - 
[T5] . describing a magnetic system coupled to an STM 
tip (tip) and a substrate (sub). The scattering region 
containing the magnetic nanostructure consists of a mag- 
netic dimer, where each atom A carries a quantum me- 
chanical spin Sa and it is characterized by an on-site 
energy eg- We assume that the tip and substrate can 
only couple to one atom at a time in the scattering re- 
gion thus to broaden the electronic level Eq by Ptip/sub 
(P — Ptip + Psub is the total broadening). We model 
the spin-spin interaction between the localized spins by 
a nearest neighbour Heisenberg Hamiltonian with cou- 
pling strength Jdd- Furthermore, we include interac- 
tion with an external magnetic field B. The electron- 
spin interaction Hamiltonian is constructed within the 
s-d model [20], where the transport electrons are lo- 
cally exchange-coupled to quantum spins through the 
exchange parameter Jgd- Thus, the full Hamiltonian 



for the magnetic nanostructure comprises three terms 
respectively describing the tight-binding electronic part 
(iJe), the spin part (i?sp) and the electron-spin interac- 
tion (ilo-sp): 
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The electron ladder operators cl^/cia create/annihilate 
an electron at site i with spin a {a —'[, I), //b is the Bohr 
magneton and g the gyromagnetic ratio. In equation ([3]) 
(T is a vector of Pauli matrices. 

The electron transport problem is solved by using the 
non-equilibrium Green's function scheme with a pertur- 
bative treatment of the electron-spin interaction, and 
where the perturbation parameter is a = Jsd/ P {p is the 
electron density of states at the spin sites, see Refs. [I6j- 
[18] for details). The diagonalisation of equation ^ gives 
the eigenenergies, e„, and eigenstates, \n), of the spin 
sytem in the interacting region. Then, at the second or- 
der in the perturbation expansion, the transition rates, 
Wni, between two eigenstates \n) and \l) induced by a 
non-spin-polarized current are given by 
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where the spin matrix elements S"^; = {n\S^\l) with i = 
{x, y, z} are that of a single spin in the chain that is 
coupled to the tip (we drop the index A). Furthermore 
C,{x) — x/{\ — er^l^^^) and /i,, is the chemical potential 
in 77-th lead, r\ = {tip, sub}. We also assume that the 
onsite energy is large enough for the density of states 
of the sample to remain constant in the small energy 
window of interest and therefore p = r/(eg -I- F^). The 
above transition rates can be used to evaluate the bias- 
dependent non-equilibrium population, P„, of the spin 
states, |n), as the steady state solution of the following 
equation 
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where T is the temperature and fce the Boltzmann con- 
stant. As P„ enters in the definition of the Green's func- 
tion for the transport electrons, Eq. ^k needs to be iter- 
ated self-consistently together with that of the electron 
propagator (see Ref. [TB] for details). 

At the second order the normalized current, /, flowing 
through a single spin can be expressed |16j as a function 



of the potential bias, V ^ as 

i(y) = ^ + ^ ^ \si,\^Puiy)i-(y - nr, 
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where the normalisation constant is given by the elastic 
conductance Go = ^ C *'r °"'' ) and where we have de- 
fined the quantity z_ {V - fl^k) = CC^^ - ^mk) - ({-V - 
fljnk)- In addition to the elastic contribution, the total 
normalized conductance, G{V) = 1 + Gi{V) + G2{V), 
includes two parts, respectively 
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At equilibrium P^ does not depend on V and Gi becomes 
the only contribution to G{V). This provides a conduc- 
tance step whenever the voltage coincides with inelastic 
energy transition flmk = \£m — £k\/e, being e the elec- 
tron charge. The conductance step height over the elastic 
conductance is then governed by the ratio pJ^^/T. How- 
ever in the non-equilibrium case Pk{V) is not constant 
and both Gi and G2 contribute, resulting in more com- 
plicated non-linear conductance profile. 

We now consider a model molecule undergoing the 
ESCE. This is described by a dimer of 5=1/2 spins ex- 
changed coupled through the bias-dependent exchange 
parameter JddiV) — Jdd + ^^^ [IH]- The parameter J^^ 
is the exchange in absence of an electric field, while b de- 
fines the value of the critical voltage, Vc — ±^y—J^^/b, 
at which the magnetic coupling switches from ferromag- 
netic (Jdd > 0) to anti-ferromagnetic (Jdd < 0). Here 
we consider, in line with density functional theory pre- 
dictions Unj, J^d = 6 meV and 6 = -0.04 meV/(mV)2, 
which give us Vc ~ 12.25 meV. The other parameters 
of the model are fixed either from experimental data or 
from theoretical considerations for systems presenting a 
similar STM setup. Firstly we consider an s-d exchange 
interaction Jgd = 500 meV, as found from theory [5T] . We 
further assume that the electronic energies of the dimer 
atoms are of the order of 1 eV, far enough from the Fermi 
energy. This produces a constant density of states p in 
the low-energy window of interest to the transport. The 
atoms are coupled to the underlying substrate strongly 
enough to produce a broadening of the onsite energies, 
Fsub = 250 meV. Furthermore, we assume that the tip is 
non spin-polarized and movable such that conductance 
spectra can be investigated for a range of Ftip. We will 
also examine the effect of magnetic fields on the spectra 
and fix the Lande g-factor to 2. All spectra are obtained 
at a temperature of 1.5 K. 

The top panel of figure [T] shows the total conductance 
calculated when the tip is positioned above one of the 



atoms forming the dimer, and the tip-to-atom coupHng is 
weak, Ftip = 5 meV. In the same panel we show how the 
spectrum is decomposed into its various parts, Gi and G2 
[see equations ^ and (IS])]. The conductance spectrum 
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FIG. 1: (Color online) Conductance spectrum for a magnetic 
dimer exchanged coupled through a voltage-dependent exchange 
parameter Jdd. (a) total normalised conductance and the different 
contributions to the total conductance. The conductance spec- 
trum is calculated by assuming a week tip-to-atom electronic cou- 
pling, Ftip = 5 meV. Note the large dip in the conductance at 
Vc ~ 12 mV, corresponding to the spin crossover transition, (b) 
Population of the singlet (black line) and triplet (red line) states 
for weak, Fjip = 5 meV (solid line), and strong, Ftip = 100 meV 
(dashed line), tip-to-atom coupling. 

displays a first conductance step at about 5 mV, corre- 
sponding to the singlet to triplet spin-excitation. Then 
at the spin crossover critical voltage, Vc, there is a second 
step in the opposite direction, i.e. the conductance dras- 
tically drops. This feature can be easily understood by 
looking at the population of the spin states as a function 
of bias, also plotted in Fig. Ill panel (b). 

For V <Vq the molecule is in the singlet state, which 
is occupied with P — \. Note that the inelastic transition 
at y ~ 5 mV does not alter the spin population since the 
current is small and not spin-polarized, so that it does not 
produce spin pumping. In other words the spin system 
always relaxes back to the ground state in between two 
inelastic scattering events. At the crossover voltage a 
change in the ground state occurs and the population of 
the singlet state rapidly decreases to zero, while that of 
the triplet states reaches 1/3. Note that the triplet state 
is three times degenerate, meaning that each of the third 
component of the S = \ state are equally occupied (there 
is no magnetic anisotropy in the model). Finally, as the 
bias is further increased a third conductance step appears 
at about 30 meV. This is now related to an inelastic 
transition between the new triplet ground state and the 
first singlet excited state, i.e. between the same spin 
states responsible for the transition found at 5 mV, whose 



energy order is now reversed. Again there is no spin 
pumping and the spin population remains unaffected. 

From the figure it is rather clear that most of the spec- 
tral features originate from Gi, which does not depend 
on the derivative of the spin population. The effect of 
G2 is evident as an additional highly non-linear contri- 
bution to the conductance at the spin crossover voltage, 
brought about by the derivative of the population of the 
spin states. 

We now investigate how the line-shape changes as the 
STM tip is brought into closer contact to the molecule, 
and Fig. [2] reports G(y) for values of Ftip in the 5- 
250 meV range. From the figure one can clearly see that 
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FIG. 2: (Color online) Conductance spectra for a spin 1/2 dimer 
undergoing the ESCE. In the upper panel we show the normalised 
conductance spectra for three different values of Ftip, while the 
contour plot in the lower panel show how the spectrum changes 
continuously between and 250 mV. Note the weakening of the 
conductance drop at Vq as the tip-to-atom coupling gets larger. 



the conductance profile gets drastically modified, in par- 
ticular in the spectral region near the spin crossover tran- 
sition. Again the results can be understood by looking 
at the evolution of the spin population as a function of 
bias, plotted for Ftip — 100 meV in Fig. [itb). The most 
striking feature is that there is a significant triplet (sin- 
glet) population for y < Vc (V < Vc), i.e. that the more 
intense current is now able to maintain the molecule in 
a steady state in which both the ground state and the 
excited state are partially occupied. It is also not sur- 
prising that the mixing is maximized at around Vc since 
the exchange coupling is reduced and the two spin states 



are less energetically separated. This feature results in 
an enhancement of the first conductance step followed by 
the associated conductance decay [SI [18] . 
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FIG. 3: (Color online) Conductance spectra for a spin 1/2 dimer 
undergoing the ESCE and subject to a magnetic field of 20 T. In 
the upper panel we show the normalised conductance spectra for 
three different values of Ftip, while the contour plot in the lower 
panel shows how the spectrum changes continuously between and 
250 mV. Note that the conductance steps corresponding to two 
inelastic transitions at both side of the spin crossover transition 
split into three equally spaced and equally intense steps. This is 
due to the Zeeman splitting of the triplet state. 

Finally we investigate the effects of applying a mag- 
netic field and for the purpose of illustration we chose 
the rather large value of 20 T (at the limit of what ex- 
perimentally achievable). Fig. [s] shows the evolution of 
the conductance spectrum as a function of Ftip. We no- 
tice that the onset of the magnetic field produces three 
equal sized and equally spaced conductance steps at 
around 5 mV. These originate from transitions between 
the 5' = ground state and the three Zeeman-split com- 
ponents of the 5 = 1 excited state. Also in this case the 
spin crossover transition produces a drastic drop in the 
conductance, which becomes less prominent as the tip- 
to-atom interaction strength gets larger. Intriguingly the 
same splitting occurs for the second conductance step at 
~30 mV, although in this case the Zeeman-split triplet 
states are resolved only for large tip-to-atom interaction. 
Note that, although for the inelastic transition at 5 mV 
the initial state is S* = and the final is S* = 1 while 
for the one at 30 mV the order is reversed, in both cases 
the associated conductance step splits with an identical 



multiplicity. This is because the ESCE simply inverts 
the order of the two spin states involved in the transition 
but does not change their symmetry. Thus the presence 
in the conductance spectrum of a drastic drop located in 
between two inelastic transitions, which split identically 
in a magnetic field, is the fingerprint of the ESCE and 
provides a unique way for identifying it. 

In conclusion we have investigated the possibility of de- 
tecting, by a spectroscopical STM experiment, the elec- 
trostatic spin crossover effect, whereby the ground state 
of a molecule changes between singlet to triplet as the 
potential voltage increases. The hallmark of the spin 
crossover transition is a drastic conductance drop, fol- 
lowing the first conductance step associated to the sin- 
glet to triplet transition. Importantly as the ground state 
changes from singlet to triplet one expects a second step 
corresponding to the opposite transition (from triplet to 
singlet). One then has the situation in which the spin 
crossover conductance drop is always placed in between 
two inelastic steps both associated to transitions between 
the singlet and the triplet. Such transitions can be eas- 
ily identified by applying a magnetic field as Zeeman 
splitting resolves the step into three equally spaced and 
equally sized steps. 
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